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Abstract 

B physics results from the Belle experiment are reviewed. Precise measurements 
of Cabbibo-Kobayashi-Maskawa matrix elements are made. Several decay modes 
are observed which will enable us to measure the CP angles, <j>( in channels other 
than B° — ► tpK° modes. New rare decay modes are observed in many channels. 
Particular attention is paid to the first observation of the electroweak penguin- 
mediated decay. 
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1 Introduction 



One of the main objectives of heavy flavor physics is the determination of the Cabbibo- 
Kobayashi-Maskawa (CKM) quark mixing matrix. || In the hadronic charged current of 
weak decay, the weak and flavor eigenstates are not identical. The CKM matrix describes 
the relation between these two sets of states. An irreducible phase in this matrix, first 
introduced by Kobayashi and Maskawa, gives rise to CP violation in the framework of 
the Standard Model (SM). 

The unitarity of the CKM matrix for b and d quark sectors leads to the expression 

V ud V* h + V cd V* b + V u V* b = 0. (1) 

This relation can be displayed as a triangle in the complex plane as shown in Figure [l|. 
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Figure 1: Unitarity triangle. 
The angles of the triangle are defined as@ 

A crucial test of the SM is to evaluate the consistency of all sides and angles of this 
unitarity triangle. This test may provide information on the dynamical origin of the 
quark mixing matrix. 
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Experimentally, the magnitudes of the elements of the CKM matrix can be determined 
by measuring decay rates (including the mixing rate) and production rates. The angles fa 
can be determined by measurements of CP asymmetries caused by the interference of two 
or more amplitudes with different weak phases. For example, sin 20i can be determined 
by the CP asymmetry induced by interference between mixing which carries the V t l phase 
and B° — > ipK° decay which carries the V* b phase as indicated by Eq. |2|. 

In this paper, we review the B physics results (except for the results which involve 
decay-time dependent analysis) from the Belle experiment. We report the determination 
of the magnitude of the CKM matrix element \V c b\ and the observation of new B decay 
modes which will lead to measurements of the CP angles fa. 

2 The Belle Detector 

The Belle detector is designed and constructed primarily to observe and measure CP viola- 
tion in B decays. Because of the high luminosity of KEKB^J {L pea k « 5.0 x 10 33 cm _2 s _1 ), 
Belle currently accumulates data at a rate of more than 4.5 fb _1 per month. This cor- 
responds to 5 million BB events per month, allowing precise measurements of B and 
charmed hadron properties. Most of the analyses presented in this paper are based on an 
integrated luminosity of 21.3 fb -1 . 

Due to the asymmetric energies of the colliding beams (3.5 GeV x 8 GeV), the T(4S) 
and its daughter B mesons are produced at (3 r y ~ 0.425 in the laboratory frame: the 
difference in B meson decay times may be measured using the difference in decay ver- 
tex positions. This key feature of the Belle experiment allows the measurement of CP 
violation through mixing in neutral B decays. 

The Belle detector consists of a silicon vertex detector (SVD), a central drift chamber 
(CDC), an aerogel Cerenkov counter (ACC), a time of flight (TOF) and trigger scin- 
tillation counter (TSC) system, an electromagnetic calorimeter (ECL), and a f^/muon 
detector (KLM). 

The SVD measures the precise position of decay vertices. It consists of three layers 
of double-sided silicon strip detectors (DSSD) in a barrel-only design and covers 86% 
of solid angle. The three layers are at radii of 3.0, 4.5 and 6 cm respectively. Impact 
parameter resolutions are measured as functions of momentum p (GeV/c) to be a\ = 
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19 2 + [50/ (p/? sin 3 / 2 6)] 2 /iin 2 and a 2 = 36 2 + [42/(p/3sin 5 / 2 6)] 2 /iin 2 , where 9 is the polar 



angle with respect to the beam direction. 

Charged tracks are primarily recognized by the CDC. The CDC covers 92% of solid 
angle in the center of mass (CM) frame, and consists of 50 cylindrical layers of drift 
cells organized into 11 super-layers each containing between three and six layers. He- 
C 2 H 6 (50/50%) gas is used to minimize multiple-Coulomb scattering. The magnetic field 
of 1.5 Tesla is chosen to optimize momentum resolution without sacrificing efficiency 
for low momentum tracks. The transverse momentum resolution for charged tracks is 
(a PT /p T ) 2 = (0.0019p T ) 2 + (0.0030//3) 2 , where p T is in GeV/c. 

Particle identification is accomplished by a combination of the ACC, the TOF and the 
CDC. The combination of these particle identification devices is a key feature of the Belle 
detector. The combined response of the three systems provide K ± identification with an 
efficiency of about 85% and a charged pion fake rate of about 10% for all momenta up to 



The CDC and ECL are used to identify electrons. The ECL also detects photons 
and measures their energy. The ECL consists of 30 cm (16.1X ) long Cesium Iodide 
(Tl) crystals. The photon energy resolution is [a E jE) 2 = (0.013) 2 + (0.0007/£) 2 + 
(0.008/E 1 / 4 ) 2 , where E is in GeV. 

The KLM is designed to detect Kls and muons. The KLM consists of 14 or 15 
modules which contain of 47 mm thick iron plates and 44 mm thick slots instrumented 
with resistive plate counters (RPC). 

In the physics analyses, we take advantage of the fact that the energy of each B 
meson is precisely known from accelerator parameters. The following kinematic variables 
are commonly used to distinguish B decay signals from backgrounds, 



where E^ e&m is the beam energy in the CM frame, p* and E* are the energies and momen- 
tum vectors of the B candidate decay daughters in the CM frame. The Mb c resolution 
is dominated by the resolution in E^ caxa and is very narrow, typically 3 MeV/c 2 . This 
variable is useful to distinguish combinatorial backgrounds because of the good resolution. 
The AE variable peaks at zero for signal events and has a typical resolution of about 10 



3.5 GeV/c. 




AE = ££*-K 



(3) 
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MeV. This variable is useful to distinguish correlated backgrounds such as feed across and 
particle misidentification since missing or extra particles, or misidentified particles cause 
a shift in the AE distribution. The M\, c distribution still peaks in the signal region for 
these backgrounds. 

In Monte Carlo simulation, the physics properties of an event are generated by the 
QQ event generator developed by the CLEO group. The detector response is simulated 
using the CERN GEANT package. [| 

A detailed description of the Belle detector can be found elsewhere. 

3 Measurement of Magnitude of CKM Matrix Ele- 
ments 

The magnitude of the CKM matrix elements, \V c b\ and \ V u b\ can be measured by studying 
semileptonic B decays. Semileptonic decay is theoretically easier to understand since there 
is no final state interactions between the lepton and hadron systems. In the naive spectator 
model, the partial decay width for the inclusive semileptonic B decay, B — > X£~V, || can 
be expressed as: 

T(B^Xrv) = ^f( 7c |K fe r + 7 u |K fe | 2 ), (4) 

where t~ — e~, \i~ and 7 9 incorporates phase space and possible strong interaction effects. 

Experimentally, the inclusive semileptonic branching fraction can be measured using 
the dilepton method introduced by ARGUS. [0 In this method, high momentum lepton is 
required to tag the flavor of one B meson and the other lepton (only electron is used) is 
used for the measurement. Primary and secondary electrons are identified by the charge 
correlation between the electron and the tagged lepton. The primary electron and the 
tagged lepton have opposite charge unless the two B mesons have the same flavor due to 
B°-B° mixing. The major backgrounds include secondary electrons from the decay chain 
b — ► c —>■ yi~V (y = s, d) of the tagged B meson and electrons from continuum e + e~ — > qq 
events, which are suppressed by using angular correlations between the tag lepton and 
the electron in the CM frame. 

Figure |2| shows the spectrum for the primary electron, dB(b — > x£~v)/dp, after sub- 
tracting the background and taking into account the effect of mixing. The branching frac- 
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Figure 2: Primary electron momentum spectrum. The histogram represents the fit result 
to the ISGW2 model. 



Table 1: V c b determination through the inclusive semileptonic branching fraction. 



Model 


\V cb \ (xlO" 2 ) 


ACCMM 


M 


4.1 ±0.1 ±0.4 


ISGWg 


I 


4.0 ±0.1 ±4 


Shifman et al. 


10| 


4.04 ±0.10 ±0.20 


Ball et al. 


n 




3.95 ±0.09 ±0.19 
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tion is obtained by a fit to the momentum spectrum predicted by theoretical models. || U 
Model dependence is very small since the measurement covers a large portion of the 
spectrum. Using a 5.1 fb _1 data sample, the inclusive semileptonic branching fraction is 
measured to be 

B(B -> Xe~V) = (10.86 ± 0.14 ± 0.47)%, (5) 

where the first error is statistical and the second error is systematic. The systematic error 
is dominated by the uncertainty in the electron identification and kinematic selection 
efficiencies. The |V^,| value is derived from the measured branching fraction and the 
model calculation of 7 C ignoring the \V u b\ term. Model dependence is relatively large as 
summarized in Table |l|. 

The value of \V c b\ can be determined with less model dependence using HQET (Heavy 
Quark Effective Theory) . [|i~2f In HQET, the partial B° — > D* + £~V decay rate is expressed 
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dTiB ° ~* d y* +l ~ V) = ^S M2 n^ M B0 ~ M D , + f\V^g{y)F{yf, (6) 

where y = v^ • vd*+ = (M| + M|,, + — q 2 )/ (2M-g M D *+ ) , v is the four-momentum divided 
by the particle mass and q 2 is the square of the four momentum transfer. The form factor 
at y = 1 (zero recoil) can be calculated with small theoretical error. Experimentally, it 
can be extrapolated from the y distribution using the following parameterization; [13[ 



g(y)F(y) 2 = j y* - l{y + if A^yf R{y) 
-2yr + r 2 2 y 

- — — (i + Ri( y y- , 

1 — ry y + 1 1 — r 



R(y) = 2 1 , %r ± T \ l + ^(y) 2 ^4) + {! + (!- R 2 (y)) y -^} 2 , (7) 



where r = M D *+/Mg . We have used a dispersion relation|14| to constrain the shape of 
the form factor, 

R^y) w i?i(l) -0.12x + 0.05x 2 , 

R 2 (y) w R 2 (l) -0.11x + 0.06x 2 , (8) 
A x {y) w A 1 (l)[l-8p 2 z+(53p 2 -15)^ 2 -(231p 2 -91)^ 3 ], 

where x = y - 1 and 2 = (y/y + 1 - Vfy/iVV + 1 + V^)- 
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Figure 3: The plot on the left shows the B° — > D* + £~V raw yield as a function of the 
measured y. The solid circles and the histogram represent the data and the fit result, 
respectively. The plot on the right shows \V cb \F(y) where the data points are derived from 
the yield by correcting for efficiency, smearing and all the factors in the differential decay 
rate. The curve displays the fit result. 



Candidate B° — > D* + l~V decays are selected by applying kinematic constraints on 
events with a electron and a D* + — > D tt + D° — > K~n + decay chain. The values of 
| V^fo | -P ( 1 ) and p 2 are extracted by a binned minimum \ 2 fit to the y distribution after 
background subtraction. The results were obtained from a 10.8 fb _1 data sample while 
fixing -Ri(l) = 1.27 and 7? 2 (1) — 0.8.|T3j Figure ^ shows the y distribution with the fit 
results. The fit yields 



V cb \F(l) = (3.54 ± 0.19 ± 0.19) x 10" 2 
p 2 = 1.35 ±0.17 ±0.18. 



Using F(l) = 0.913 ± 0.042, JTJ \V cb \ is determined to be 

= (3.88 ± 0.21 ± 0.21 ± 0.19) x 10" 2 . 



(9) 



(10) 



where the third error is theoretical. The systematic error is dominated by the uncertainty 



in the tracking efficiency for slow pions from D 



*+ 



7T 



decay. 



The magnitude of V ub is one of the smallest and poorly measured parameters of the 
CKM matrix. This is the key element needed to evaluate the consistency of the SM 
with the large CP angle, sin20i = 0.99 ± 0.14 ± 0.06, [|T6| measured by Belle. Exclusive 
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Figure 4: M^ c distribution for B° — > n + l v decay. The unshaded histogram shows the 
signal component while the hatched histograms show various background components. 



semileptonic b — > ut~v decay is one of the most promising modes to determine \V U \\ 
provided that the form factors are reliably calculated by theoretical models. 

The branching fractions for B° — > p£~T7 and B° — > n£~V have previously been mea- 
sured by CLEO using "neutrino reconstruction" technique. ||17|| In this method, the four 
momentum of the undetected neutrino is inferred from the missing momentum and en- 
ergy in an event. The missing momentum is a poor measure of the neutrino momentum, 
when the event has multiple neutrinos, missing or extra particles. Such events are re- 
duced by requirements on total charge and the number of leptons. The missing mass- 
squared of the event in the CM frame, (M miss ) 2 = (E^ iss ) 2 — (p^iss) 2 ; i s required to satisfy 
|(M miss ) 2 | < 2 (GeV/c 2 ) 2 to further suppress events with poorly reconstructed neutri- 
nos. The major background from b — > c£~T> is reduced by requiring \p\\ > 1.2 GeV/c, 
\p}\ + \p*\ > 3.1 GeV/c and other kinematic constraints. We require |Ai?| < 0.3 GeV 
to reduce feed-across as well as combinatorial backgrounds. The signal yield is ex- 
tracted by a fit to the M^ c distribution. A Monte Carlo (MC) simulation is used to 
provide both signal and background shapes. We find 107 ± 16 B° — > tt + £~V events in 
a 21.3 fb _1 data sample as shown in Figure f|, corresponding to a branching fraction of 
B(B° -> n + £-v) = (1.24 ±0.20 ±0.26) x 10~ 4 . The systematic error due to uncertainty of 
the neutrino finding is dominant. Further studies are required to extract \V u b\ from this 
measurement. 
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Figure 5: M bc distribution for 5° — > £)* ± D =F and 5° — ► D* + D*~ modes. The solid line 
shows the sum of signal and background while the dashed line shows background. 



4 Toward 6, Measurements 



The main motivation of the Belle experiment is to measure all the CP angles (<pi). Such 
measurements are essential to probe new physics and overconstrain the CKM matrix. We 
report some measurements and observations of new decay modes which demonstrate our 
capabilities to measure these CP angles. 

Doubly-charmed decay modes B° — > D* + D~ , D + D*~ (referred as D^D^ hereafter), 
D*+D*- can foe used to measure^] <fii in addition to the gold plated mode B° — > ipK°. 
The CP asymmetry parameter measured in this mode could deviate from the expected 
value due to a sizable penguin contribution, which may provide evidence for new physics. 
The B° — > D* ± D T decay mode has not been observed so far. 

Candidates are selected using the decay chains, D~ — > K + tt~tt~ and D* + — > D°7r + , 
D° — > K~ir + , K~7i + 7r°, K~ tx+ tx~ . We also use the D° — > KgiT + TT^ , Kg — > 7r + 7r~ decay 
chain for the B° — > D* + D*~ mode. The signal yield is extracted by a maximum likelihood 
(ML) fit to the Mb c distribution for events within the AE signal region (|Ai£| < 20 MeV). 



A Gaussian and an ARGUS functions^] are used to represent the signal and background 
shapes. The fits yield 11.2±4.0 B° -> D*^D T events and ll.0i3.7B -> D* + D*~ events 
in a 21.3 ftr 1 data sample with statistical significances of 4.1<r and 5.0a, respectively. 
The statistical significance is defined as ^—2 ln(£ /£ max ), where £ max is the maximum 
likelihood and £q is the likelihood values at zero signal yield. Figure |5] shows the M^ c 
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Figure 6: cos a distribution for partially reconstructed B° — ► D* ± /} =F decay. 

distributions with the fit results. The branching fractions are measured to be 

B(B° — > D* ± D =F ) = (1.04 ± 0.38 ± 0.22) x 10" 3 , 
B(B°^D* + D*-) = (1.21 ± 0.41 ± 0.27) x 10" 3 . (11) 

The systematic error is dominated by the uncertainties in tracking efficiency and fit 
method dependence. 

The B° — > D* ± Z} =F decay mode is also observed using a D* + partial reconstruction 
technique. In this method, after the D* + — > D°7r + decay we do not reconstruct the 
subsequent D° decay in order to increase the overall detection efficiency. Charmed mesons 
from B° — > D* ± D =F are almost back-to-back in the CM frame since 5 mesons are produced 
almost at rest. A slow pion from the D* + — ► D°tt + decay approximately retains the 
momentum direction of the parent D* + because of the small energy release in this decay. 
Thus the angle a between the slow pion and the D~ are almost back-to-back and can 
be employed as a signature. However, the partial reconstruction method suffers from a 
relatively large background. In particular, continuum e + e~ —> cc events can produce D* + 
and D~ in a back-to-back configuration, which results in similar angular correlation. The 
partial reconstruction method also introduces a difficulty for the <f>\ measurement since the 
particles which are not used for the reconstruction complicate the vertex reconstruction 
and the flavor tagging of the accompanying B meson. These problems can be alleviated 
by requiring a high momentum lepton (p* t > 1.1 GeV/c) in the event. This requirement 
heavily suppress the continuum background, and at the same time, provides a clean flavor 
tag and a good vertex reconstruction for the <pi measurement. The remaining backgrounds 
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can be further suppressed by exploiting the angular correlation between the tag lepton 
and D~. Figure |6| shows the cos a distribution observed in a 21.3 fb" 1 data sample. 

A fit to the cos a distribution for the lepton tagged sample yields 35.8 ± 11.3 signal 
events. The signal and background shapes are obtained from the MC simulation except 
for the background from Cabbibo-favored B° — > D* + D~ decays for which the data is used 
to estimate the amount and the shape of the background. The branching fractions are 
measured to be 

B(B° -> D^D T ) = (1.78 ± 0.56j£re) x 10 ~ 3 - ( 12 ) 

Uncertainties in the background shape are the dominant source of the systematic error. 
This measurement clearly demonstrates that the partial reconstruction technique can be 
used for a </>i measurement with B° — ► _D* ± .D =F and B° — ► D* + D*~ decays. 

The gluonic penguin decay mode B° — > <fiK s can also provide an independent mea- 



surement of the CP angle, 0i.|2"O|] This decay mode proceeds through a b — > sss transition, 
which is forbidden at the tree level in the SM, but are induced by second order loop dia- 
grams (penguin or box diagrams). The CP asymmetry parameter measured in this mode 
is of special interest since it is sensitive to the possible exchange of non-SM particles in 
the loop|2]J and may deviate from the expected value. 

The B° — > <pK s candidates are reconstructed via — > K + K~ and K s — > it + tt~ 
decays. The dominant background arises from continuum events. Event shape variables 
are used to suppress the continuum background. The most powerful variable is the cosine 
of the angle between the B candidate thrust axis and the thrust axis of the rest the 
event (cos#t)- The cos^t distribution for the signal is flat while it is peaked at ±1 for 
continuum events. This variable is combined with other variables such as the B flight 
direction and the helicity angle for pseudoscalar- vector final states {4>K S , <pK~) using a 
likelihood ratio C1Z = £s/(£s + £b) where Cs an d C-b denote the signal and background 
likelihoods. The likelihood is the product of probability density function in each of the 
discriminating variables. 

An extended unbinned ML fit is performed in AE and M^ c simultaneously to extract 
the signal yield. In the extended ML fit, the sum of the signal and background yield is 
allowed to be different from the total number of the event in the fit. The signal shape 
is represented by Gaussian function in both AE and M^ c . The background shape is 
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Figure 7: Mb c and AE distributions for B° — > 4>K S decay. 



represented by a linear function in AE and an ARGUS function in M^ c . The fit yields 
8. Oilg B° — > 4>Kg events with a statistical significance of 4.2a in a 21.3 fb _1 data sample. 
Figure ^ shows the M\ )C and AE distributions with the projection of the fit result. We 
have also observed other B — > <pK modes and obtained the branching fractions as 



(0.89^ ±0.10) x 10" 5 , 
(1.30^ ±0.21) x 10" 5 , 



(13) 



'i i O+0.22 
k ±.±Z_q 20 



±0.14) X 10' 



Here, K*° refers to the K*°(892). The dominant systematic errors come from uncertainties 
due to tracking and Kg reconstruction. 

Charmless B decay B° — > 7t~7t + is one of the most promising modes to measure the 



CP angle 2 . We also need to study B c 



n°7i° and B' 



7T 7T 



decays to disentangle 
the effect of penguins. |E| Experimentally, K/ir separation at high momentum is the 
key to distinguish the nir mode from the Kit mode. The ACC detector in the Belle 
plays an essential role in this regard. The K/tc separation capability is measured using 
kinematicaly selected D* + — > D°tt + , D° — > K~n + decays in the data. The efficiency 
is measured to be 92% for pion and 85% for kaon. The misidentification probability is 
measured to be 4% for pion (true pion fakes kaon) and 10% for kaon. The dominant 
background comes from continuum events. The Super-Fox- Wolfram (SFW) variable 



which is an extension of the normalized Fox- Wolfram moment s [g4j] is employed to suppress 
the continuum background in addition to the cos#t- These event shape variables are 
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Figure 8: AE distribution for B° — > tx tt + decay. The fit function and its signal and 
cross-talk components are shown by the solid, dashed and dotted curve, respectively. 

combined into a Fisher discriminant^] since they are correlated. The resulting Fisher 
discriminant is combined with the B flight direction using a likelihood ratio. This selection 
rejects 95% of the continuum background while retaining 40% to 50% of the signal. 

The signal yield for B° — > 7t~tt + mode is determined from a fit to the AE distribution 
where there is kinematic separation between tt~tc + and K~n + modes. The signal shape 
is modeled by a Gaussian function while the background shape is modeled by a linear 
function and a Gaussian function. The Gaussian function is introduced to account for 
the background from Kir mode. The mean value of this background is shifted by about 
—50 MeV since the kaon is misidentified as pion. The normalization of all components 
are free parameters in the fits. The Gaussian background yield extracted from the fit 
is consistent with the misidentification probability mentioned above. Figure |8| shows 
the AE distribution for B° — > 7r~7r + mode along with the fit result. The signal yield for 
B~ — > 7r~7T° mode is determined from the simultaneous fit to the AE and Mb c distribution 
because of the long tail for the AE signal shape. We find 17.7^4 events for the B° — > 
7Y~7i + mode and 10. 41^3 events for the B~ — > n^n mode in 10.4 fb" 1 of data. The 
branching fractions are measured to be 

B{W^W+) = (0.56±gS±gS) x 10- 5 , 

B(B- - TT-vr ) = (0.78^;?2) x 10- 5 . (14) 



The interference between B — > D°K {b cus) and B — > D°K (b — > ucs) 
decays provides a theoretically clean determination of the angle 03. |26|, 27] When D° and 
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D° decays into a common CP eigenstate mode, the CP angle 03 can be related to the 



following observables||26|| assuming no D°-D° mixing; 



A EE 



R 



T(B- 



DcpK- 



T(B- 



DcpK- 



T(B~ -> D CP K-) + -> £> CP tf+) 

2£yr sin 5 sin 3 , 

p' / p = 1 + r 2 + 2£/r cos 5 cos 03, 

r(£r - p C p^) + r(p+ - p C p^ + ) 



(15) 



P = 



r(B~ — > Dcpi?-) + r(£?+ 
r(p- -> pAfc-) + r(p+ - 



-> D(jp7T' 

D°K+) 



D 



T(B- 



T(B- 

A(B~ ^D°K-)/A{B 



-> P 7T 



where Pcp is the CP eigenstate of the _D° meson, 5 is the strong phase-difference between 
B~ -> Z^A"" and S~ -> P>°A-, and is the CP eigenvalue of D CP . 

We have studied P~ — > D°K~ decays where D° decays into K~ir + or into a CP = +1 
eigenstate (K~K + , 7r~7r + ). Continuum events are the dominant background and are re- 
duced by a selection based on event shape variables. Tight particle identification is applied 
to reduce the background from B~ — > D°tt~ decays. The signal yields are extracted from 
a fit to the AE distribution that accounts for the remaining B~ — > P°7r~ background. 
The background function includes an ARGUS function to model the combinatorial back- 
ground and a Gaussian function to model the B~ — > D°tt~ background. Figure [| shows 
the AE distributions with the fit results. It should be noted that peak position for the 
signal is shifted by —49 MeV since AE is calculated with the assumption of a pion mass 



for the prompt kaon. We find 1278 ± 37, 114 ± 11 and 70 ± 11 B' 



D°ii events 



and 85 ± 10, 12.3 ± 3.9 and 4.9 ± 5.4 B~ -> D°K~ events with subsequent D° decays 
into K~ir + , K~~K + and 7r~7r + modes, respectively, in 21.3 fb _1 of data. The statistical 
significance of the B~ — ► D°K~, D° — > K~K + signal is 4.3cr. These results gives 



.4 
R 



0.04^° ±0.15, 
1.39 ±0.53 ±0.26. 



(16) 



This measurement is an important first step toward a 03 measurement. 
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Figure 9: AE distribution for for B — ► decay with subsequent Z) decay into 

K~7i + , K~K + and 7r~7r + . The fit function is described in the text. 

5 Charmless Hadronic Decays 



Charmless hadronic B decays proceed primarily through b — > u tree diagrams and b — > s 
penguin diagrams, which provide a rich ground to search for direct CP violation, although 
it is theoretically difficult to relate the CP angles to the CP symmetries measured in 
these decay modes. In addition, there are possibilities to extract CP angles from ratios of 
B — > Kit decay rates. 0] Furthermore, the involvement of the penguin diagrams makes 
these decays sensitive to the effect of new particles at higher mass scales in the loop. 



It is a theoretical challenge to explain the unexpectedly large branching fractions for 
B — > rj'K and B —>■ rjK* decays reported by CLEO.[^] Since it may suggest new physics 
beyond the SM, it is necessary to experimentally confirm the CLEO results and search for 
direct CP violation which may be enhanced by new physics. The large branching fraction 
for B° — > rj'Kg also opens the possibility of measuring mixing- induced CP violation. 

We analyze B~ — > f]'K~, B° — > rj'Kg and B° — > rjK* decay through the decay 
chains 7/ — > (77 — > 77), rf — > p°7 (p° — > 7T + 7r~), Kg — > 7r + 7r _ and — > fT _ 7r + . 

We also use the 77 — > 7r + 7r~7r° decay for the 5° — > 7/i^* analysis. After background 
suppression using variables such as cos^t, SFW variable, B flight direction and helicity 
angle for B° — > r/i^* mode, we perform extended ML fits to both Mt, c and AE to extract 
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Figure 10: M^ c and AE distributions for B° — > r/i^s mode. Solid and dashed curves 
show the projections of the signal+background and background functions. 



signal yields. We observe 71.4 B~ — > rj'K" and 16.5 B° — > rfKg events with statistical 
significances of 12. Oct and 5.4<r, respectively, in 10.1 fb _1 of data, and also 22.1 B° — > i]K*° 
events at 5.1a in 21.3 fb _1 of data. Figure [H] shows the M^ c and AE distributions for 
B° — > rj'Kg mode with the projections of the fit function. The branching fractions are 
measured to be 



B(B~ -> ri'K-) 
B(B° -> r/Z°) 
£(5° -> ryiT ) 



(79ig ± 9) x 10~ 6 , 
(55±i| ± 8) x 10 

j+5.4 



-6 



(21.2+^ ±2.0) x 10" 



(17) 



These results confirms the large branching fractions in these modes observed by CLEO. 
In addition, we have measured the decay asymmetry for the B~ — > i]'K~ mode, 

T(B~ -> rj'K~) - T(B + -> 77'^+) 



.4 



CP 



F(B~ -> r/'K-) + r(£+ 
+0.06 ±0.15 ±0.01. 



The result is consistent with no direct CP violation and with the SM prediction. 

Three-body charmless decays B~ — > K~h + h~ (h refers to a charged kaon or pion) 
proceed through variety of diagrams such as b — > uud, b — > ccs and b — > sss, which 
provide an excellent environment to search for CP violation due the interference of these 
amplitudes. In particular, a great deal of attention is paid to the interference between the 
former two amplitudes because it is sensitive to the CP angle 03. The B~ — > 7t~7t + tt~ 
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Figure 11: AE distributions for B~ — > K~7r + 7r~ and B~ —> K~K + K~ modes. Points 
are data, histograms are backgrounds estimated by the MC, and solid and dashed curves 
show the signal+background and background functions. 



decay is a good example since the non-resonant mode can interfere with the decay B — > 

In this analysis, we reconstruct the decays B~ — > K~h + h~ without any assumption 
on the intermediate hadronic resonance. As in other rare B decay modes, continuum 
events are the dominant background source and are suppressed using various event shape 
and kinematic variables. In the B~ — > K~tt + tt~ mode, we have large backgrounds from 
Br -> D°n-, D° -> ^-tt+ and Br -> tf-^, V C) -> A* + A*~ 0° refers to J/V> and 
t/>(2S)) where muons are misidentified as pions. These backgrounds are suppressed by 
requiring \M K - n + - M D o\ > 0.1 GeV/c 2 , \M h+h - - Mj^\ > 0.07 GeV/c 2 and \M h+h - - 
^V(2S)| > 0.05 GeV/c 2 . In the B~ — > K~K + K~ mode, the background from the — > 
D°K~, D° -> decay is rejected by requiring |M^- X + - M D o| > 0.025 GeV/c 2 . 

Signal yields are obtained from fits to the AE distributions. We find 177 ± 20 B~ — > 
K~7t + tt~ events and 162 ± 16 B~ — > K~K + K~ events in a 21.3 fb _1 data sample, which 
corresponds to branching fractions; 

B{B~ -> K-n + n-) = (58.5 ± 7.1 ± 8.8) x 10" 6 , 
B(B~ —> K~K + K~) = (37.0 ± 3.9 ± 4.4) x 10" 6 . (19) 

Figure [11] shows the AE distribution for the B~ — > i^~7r + 7r~ and the B~ — > K~ K + K~ 
decays after the background suppression. Our result is higher than the upper limit, 
B(B~ -> K-7T+7T-) < 28 x 10~ 6 , reported by CLEO.^l The CLEO's analysis required 
the mass for any pair of particles to be above 2 GeV/c 2 , which effectively eliminates the 
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Figure 12: The 7r + 7r~ mass distribution in the B~ — > K~tt + tt~ decay and the K + K~ mass 
distribution in the K~K + K~ decay. The signal function is modeled by a Breit-Wigner 
function convolved with a Gaussian. 



low mass resonances that are found to dominate the signal. 

Further studies of intermediate resonant states of these decays are made with a 
Dalitz plot style analysis. Clear contributions from B~ — > K*°tt~, K*° — > K~ir + and 
B~ — > i^~/o(980), /o ~~ > 7T + 7r _ , and no K~p° signal are observed in the B~ — > K^-k + -k^ 
decay. We also find broad resonances in K~n + and 7r + 7r~ mass around 1.4 GeV/c 2 and 
1.3 GeV/c 2 , respectively In the B~ — > K~K + K~ mode, we find a clear 0(1020) reso- 
nance and a broad K + K~ enhancement around 1.6 GeV/c 2 . Larger statistics is required 
to identify these resonances using a Dalitz plot analysis with interference terms between 
the resonances. 

We also find clear signals for B~ — > XcqK~ , Xco — > 7r + 7r~, K + K~ in B~ — > K~tt + tt~, 
K~K + K~ decays. Figure shows the 7r + 7r~ and K + K~ mass distributions in B~ — > 
K~it + -k~, K~K + K~ decays. The fits to these distributions yield 15.5jl|;| events in the 
7r + 7r~ mode and 7.7lg^ events in the K + K~ mode at statistical significances of 4.8a and 
3.2cr, respectively. The fit also gives a sizable peak shift in the K + K~ mass spectrum, 
which may be due to interferences from other K~K + K~ final states. Because of this 
uncertainty, we determine a branching fraction from the 7r + 7r~ mode only, 

B{B-^Xc K-) = (8.0±l;t ± 1.9 ± 1.1) x 10~ 6 , (20) 

where the third error is due to the uncertainty of the Xco vr + 7r _ branching fraction. 
This observation provides evidence for a significant nonfactorizable contribution in B to 
charmonium decay process and suggests the existence of a sizable B~ — > Xco 71 " - decay, 
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which could be used for a 3 measurement. 



6 Color-Suppressed Decays 

Color-suppressed B decays such as B° — > D^*'°h°, where h° refers a light neutral mesons, 
proceed through an internal spectator diagram which is suppressed by color-matching. 
Since the branching fractions for these modes are expected to be very small [32] (< 1CT 4 ), 
studies of color-suppressed decay could provide useful information on final state interac- 
tions and hadronic -B-decay models. 

In this study, B° -> D^°h° decays are reconstructed through D*° — > D°h°, D° — > 
K~7i + , K~7[ + 7t°, K~7t + 7t~7t + where h° is reconstructed through 7r° — > 77, rj — > 77, 
7r + 7r _ 7T° and — > 7r + 7r~7r° modes. In addition to the continuum background which is 
suppressed by event shape variables and helicity angles, we pay attention to the back- 
ground from color-favored decays. The B° — > D* + p~ mode can give the same final state 
as D°r) and D°u modes. Such background is mostly removed by the 77 and u mass con- 
straints and its contribution can be evaluated from the rj and u mass sidebands. The 
B~ -> D^°p- decays can contaminate the B° — > D^*)°7r mode if the momentum of the 
tc~ from p~ is very small. These backgrounds are reduced by rejecting the events that 
can be reconstructed as B~ — > D^°p~. These events also have AE values lower than the 
signal due to the missing pion. 

Signal yields are extracted by fits to the AE distributions taking into account the 
backgrounds from color-favored modes. We observe 126±16 events in D°n° mode, 26.41^ 
events in the D* n mode, 22. 1^3 events in the D°r) mode and 32.5^g;g events in the D°uj 
mode in a 21.3 fb _1 data sample corresponding to statistical significances of 9.3a, 4.1cr, 
4.2cr and 4.4c, respectively. Figure |TB| shows the AE distributions for the D tt and D°uj 
modes along with the fit results. The branching fractions are measured to be 

B(B° -> D\°) = (3.1 ±0.4 ±0.5) x 10" 4 , 

B(B°^D* n°) = (2.T±8;|±g;g) x 10" 4 , 

B(B°^D°ri) = (1.4+^ ± 0.3) x 10" 4 , (21) 

B(B° -> D°uj) = (1.8 ± 0.5±°; 4 ) x 10~ 4 . 

These results are consistently higher than the predictions, (0.5 ~ 1.0) x 10 -4 , by a fac- 
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Figure 13: AE distributions for B° — > D tt and B° — > D°u modes. The solid lines 
show the fit functions, and the dashed lines show the sum of the feed-across and the 
combinatorial background, with the latter shown separately as the dotted lines. 



torization model, [32] which could indicate the presence of non-factorizable effects such as 
final state interactions, or some corrections to factorization. 



7 Penguin- Diagram Mediated Decays 

The b — > s transition is a penguin- diagram mediated FCNC process. Since it is forbidden 
at tree level in the SM, its small amplitude make it sensitive to effects caused by exchange 
of non-SM particles in the penguin loop. 

The CLEO group reported the first observation of the B — > X s 7 radiative penguin 
decay. p3| Here, B represents B° or B~ mesons and X s represents a hadron system that 
includes a s quark. The measured branching fraction for this process provides the most 



stringent indirect limit on the charged Higgs mass range ||3~4|| and a constraint on the mag- 



nitude of the effective Wilson coefficient of the electromagnetic penguin operator [ | . 
The electromagnetic penguin decays B —>■ X s l + l~~ are essential to further constrain the 
Wilson coefficients. The magnitude and the phase of the coefficients Cy ff , C| ff , C w can be 
determined by measuring the dilepton invariant mass distributions and forward-backward 
charge asymmetry of the dilepton and the B — > X s ^ rate. ||35| These measurements are 
crucial to obtain definitive evidence for new physics. 

The branching fraction for B — > X s "f has been calculated to 10% precision in the 



framework of the SM including next-to-leading order QCD corrections. [36] It is important 



21 



to measure the branching fraction to a precision of 10% or better to explore or limit non- 
SM theories. 

The inclusive B — * X s ^ decay is reconstructed by combining the X s system with a 
photon. The X s system is formed by combining a charged kaon or Kg with 0-4 pions 
which may include one tt°. Combinatorial background is reduced by requiring Mx s < 
2.05 GeV/c 2 . The SFW variable is employed to suppress continuum background. In 
addition, the SFW sideband is used to model the background shape for the continuum 
events. The fit to the M^ c distribution yields 107 ± 17 events in a 5.8 fb _1 data sample, 
corresponding to a branching fraction of 

B(B -> X sl ) = (3.36 ± 0.53 ± 0.42l°j°) x 10" 4 . (22) 

The third error is due to extrapolation of the Mx s spectrum from the region Mx a < 
2.05 GeV/c 2 . Our result is consistent with the SM prediction |6| of (3.28 ± 0.33) x 10~ 4 . 



In order to reduce the systematic error due to the Mx a spectrum, we need to under- 
stand the resonant structure of the X s system beyond the B — > K* ^ decay which has 
been well measured. The experimental search for higher kaonic resonances yielded only 
an indication of B — > .fQ (1430)7 so f ar - 111 this analysis, we study the kaonic resonances 
which decay into K~n + , Kgn~ , K~ir°, and K^-k + -k^ modes. 

In the B° — > K~7r +/ y mode where the K~it + mass is required to be consistent with 
K*°(U30), a fit to the M bc distribution yields 29.1 ± 6.7 events in a 21.3 fb -1 data 
sample. The helicity angle distribution is analyzed to distinguish the i^2°(1430) signal 
from i^*°(1410) and non-resonant modes. A fit to the helicity angle distribution yields 
20.1 ± 10.5 events for the K 2° {14:30) component, which leads to a branching fraction of 

B(B° -> 7q°(1430)7) = (1.26 ±0.66 ±0.10) x 10" 5 . (23) 

We reconstruct B~ — > K~7T + 7T~'-f decays in B~ — > K*°7i^ / ~f and B~ — > K~p 0/ y modes 
where K~n + or 7r + 7r~ mass is required to be consistent with K*° or p°. M K - n + w - is 
required to be less than 2.0 GeV/c 2 to reduce combinatorial backgrounds. Fits to the 
Mb c distributions yield 46.4 ±7.3 events in the B~ — > i^* 7r~7 mode and 24.5 ±6.4 events 
in the B~ — > K~p 0/ ~f modes. Figure O shows the Mb c distributions for B~ — > K* 7T~^f 



and B — > K p°7 modes with the fit results. After subtracting the backgrounds such as 
B~ — > K~7T + 7T~~f non-resonant decays, feed-across and other B —>■ Xgj decays, we obtain 
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Figure 14: Mb c distributions for B — ► A'*°7r 7 and B —> K p°7 modes. The solid lines 
show the fit functions, and the dashed lines show the background. 

signal yields of 39.7 ± 7.4 events in the B~ — > A*°7r~7 mode and 22.2 ± 6.5 events in the 
B~ — > K~ p°7 modes, which gives branching fractions of 

B(B~ K*\- T , M^o^ < 2.0 GeV/c 2 ) = (5.6 ± 1.1 ± 0.9) x lO" 5 , 
B(B~ -> K-p° 7 ; M^-^o < 2.0 GeV/c 2 ) = (6.5 ± 1.7+^ x 10~ 5 . (24) 

The sum of measured exclusive modes accounts for about a half of the total B — * X s ^ 
branching fraction. 

A great deal of attention has been paid to the B — > X s £ + £~ decay since it provides 
information necessary to determine the Wilson coefficients, Cy , C| ff , and C\q. Some non- 
SM models such as SUSY gives significantly different values for these coefficients from the 
SM prediction due to non-SM particle contributions to the loop in the penguin diagram. 
However, no experimental evidence has been observed for the B — > X s £ + £~ decays. 

We have searched for B — > X s £ + £~ decay using both an exclusive and inclusive ap- 
proach. The background from B ipC^K, ipO — > £ + £~ is rejected by requiring the 
dilepton invariant mass to be outside of ipO mass region. We have a wider veto region 
for the electron mode because the electron tends to lose its energy due to initial state 
radiation or bremsstrahlung. Several event shape variables are combined into a Fisher 
discriminant to suppress continuum background. The missing energy of the event is used 
to suppress the major background from semileptonic B decays. The Fisher discriminant 
variable and the missing energy are combined with kinematic variables into likelihood 
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Figure 15: M^ c distributions for B — ► K£ + £~ modes. The solid lines show the fit functions, 
and the dashed lines show the background in the left column. The solid lines show a fit 
to the ARGUS function in the right column. 

ratios to maximize the background rejection capability. We reject 85% of continuum 
background and 45-55% of BB background while retaining 70-75% of the signal. 

In the inclusive analysis, the signal yield is obtained from a unbinned ML fit to the 
M bc distribution with the sum of a Gaussian function for the signal and an ARGUS 
function for the background. The background shape is determined by the fit while the 
signal shape is calibrated using the B — > J/ipK sample in data. We find 11.4lf;g events 
with a statistical significance of 2.7a for the mode in a 29.5 fb _1 data sample. 

In the exclusive analysis, the background shape is determined using a MC sample 
corresponding to 400 fb _1 due to lack of statistics in the data sample. We also take into 
account the background due to lepton misidentification where B — > K^h + h~ decay is 
misidentified as B -»• K { *h + r where represents K~ , K°, K*~(892) and K*°(892). 
This background is included as a small additional Gaussian background function in the 
fit. We find 9.5l|;f events in the Kfi + fx~ mode, 4.ll|;^ events in the Ke + e~ mode and 
6.3l3;o events in the K*e + e~ mode with statistical significances of 4.7a, 2.5a and 2.5a, 
respectively. We do not observe a significant signal in the K*fi + fi~ mode. The probabil- 
ity of an upward fluctuation of the background to the observed number of events in the 
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Kfi + fi mode is 5.5 x 10 6 , which corresponds to 4.4<r for a Gaussian probability distri- 
bution. When we combine the Kp^pT and Ke + e~ modes, we observe I3I3J events with 

K£ + £- 
modes with the fit 



a statistical significance of 5.5<r further establishing the observation of the B 
decay. Figure [15] shows the M^ c distributions for the B — > K£ + £ 
results. The left column shows the data and the right column shows the background 
distributions in the MC sample. The branching fraction is measured to be 



B(B 



(0.75±g;|f ± 0.09) x 10" 



(25) 



This result is consistent with some of the SM predictions |37] which cover the range (0.42 



0.57) x 10 6 , although one model|38| gives a slightly lower value of (0.33 ± 0.07) x 10~ 



8 Conclusions 

We have established the existence of the electroweak penguin-mediated decay B — > 
K£ + £~, which opens a new window to search for physics beyond the SM. We have also 
made the following new observations. 

• Color-suppressed B decays, B° — > D°n°, D 0r q, D°uj. 

• Three-boy charmless B decays, B° — > K~h + h~, in which B~ — > XcoK~ decay is 
observed. The large branching fractions for XcoK~ decay and the color-suppressed 
decays pose challenges to factorization models. 

• B~ — > DcpK~ decay. We have made an asymmetry measurement, which is the 
first step toward a 03 measurement. 

• B° — > D* ± D =F decay with full and partial reconstruction technique. This is another 
promising mode that can be used to measure <p\. 

• B~ — > K* 7T~^f and B~ — > K~p 0/ y. The sum of the exclusive branching fractions 
accounts for about 50% of the inclusive branching fraction. 

In addition to the above new observations, we have observed and measured branching 
fractions for B — > 4>K, B° — ► D* + D*~ and B — > nir decays. These modes will provide 
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measurements of 0i and 2 in the near future. We have also confirmed the large branching 
fractions for B — > rj'K and i? — > r]K* decays, which is another theoretical challenge. 

We have made precise measurements of |V^,| via inclusive semileptonic B decay and 
B° — > D* + l~V decay. We have measured a branching fraction for 5° — > 7r + £ _ z7 decay, 
which will lead to a |T4b| measurement. A precise measurement of |T4&| is the key to 
evaluate the consistency of the Standard Model with the large sin 20i value measured by 
the Belle experiment. 

In conclusion, we have demonstrated that we are ready to take the next step. We can 
measure all the CP angles in many different modes in addition to measuring sin20i via 
the B° — > ipK° mode, and we have many effective tools to probe new physics beyond the 
Standard Model. 
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